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We have constructed an optical clock with a fractional frequency inaccuracy of 8.6 x 10~ 18 , 
based on quantum logic spectroscopy of an Al + ion. A simultaneously trapped Mg + ion serves 
to sympathetically laser-cool the Al + ion and detect its quantum state. The frequency of the 
1 So'f->' 3 Po clock transition is compared to that of a previously constructed Al + optical clock with 
a statistical measurement uncertainty of 7.0 x 10~ 18 . The two clocks exhibit a relative stability of 
2.8 x 10 _15 t _1//2 , and a fractional frequency difference of —1.8 x 10~ 17 , consistent with the accuracy 
limit of the older clock. 



Optical clocks based on petahertz (10 15 Hz) transi- 
tions in isolated atoms have demonstrated significant im- 
provements over the current cesium primary-frequency- 
standards at 9.2 GHz. They also shed light on funda- 
mental physics, such as the possible variation of physical 
constants. While the merits of laser-cooled ion optical 
frequency standards were known [l], further devel- 
opments were required to permit their use. Sub-hertz 
linewidth lasers [3[ enabled single ions to be probed with 
sufficient resolution for high-stability clock operation, 
and control of external-field perturbations allowed such 
clocks to operate with an inaccuracy below 10 -16 [3, HI]. 
For comparison, cesium standards that realize the SI sec- 
ond have reached an inaccuracy of 3 x 1CP 16 [f|, and 
an optical lattice clock based on Sr atoms has been re- 
ported 0] with an inaccuracy of 1.5 x 10~ 16 . Here we de- 
scribe an Al + ion clock with an inaccuracy of 8.6 x 10 -18 . 

The 1 S o 3 P transition in A1+ at 1.121 PHz is of in- 
terest due to its low sensitivity to electromagnetic pertur- 
bations and its narrow natural linewidth of 8 mHz [H, Q . 
Al + has the smallest sensitivity to blackbody radiation 
loL 11 1 among atomic species currently under considera- 
tion for clocks, thus relaxing the requirement on ambient 
temperature control. However, the absence of an accessi- 
ble allowed optical transition prevents the internal state 
of Al + ions from being detected by conventional meth- 
ods, and the ion cannot currently be directly laser-cooled. 
Quantum logic spectroscopy (QLS) [l2| overcomes these 
difficulties by trapping a "logic ion" that can be directly 
laser-cooled together with the A1+ clock ion. The coupled 
motion of the two ions allows for sympathetic cooling, as 
well as the transfer of the clock ion's quantum state to 
the logic ion, where the state can be measured. 

The clock described here shares features with our 
previously-constructed Al-Bc clock [B], Q, but also in- 
cludes many changes, making a comparison of the two 
clocks a valuable test of systematic errors. In the new 
clock, the 9 Bc + logic ion has been replaced by 25 Mg + , 
whose mass closely matches that of 27 A1 + . Laser-cooling 
inefficiencies due to mass-mismatch are thus suppressed. 
The Al-Mg ion trap is a linear Paul trap built from 
tool-machined all-metal electrodes (Fig. 1). This con- 
struction differs from the Al-Be trap that was built from 



1070 nm 
fiber laser 




Frequency 
recording 



FIG. 1: Setup for comparing the frequencies of the two Al + 
clocks. The 4 th harmonic of a fiber laser is locked to the 
Al + 1 So<-> 3 Po clock transition in the Al-Mg apparatus with 
a fixed offset frequency (applied to AOM1). Another laser 
beam derived from the same laser probes the clock transition 
in the Al-Be apparatus, where the laser frequency is locked to 
the clock transition in a separate digital feedback loop that 
controls AOM2. The record of frequencies applied to AOM2 
represents the difference in clock frequencies. BS: beam split- 
ter; AOM: acousto-optic modulator; x2: frequency doubler. 



laser- machined and gold-coated alumina electrodes [1 31 ] . 
and it exhibits reduced RF-micromotion-inducing elec- 
tric fields. 

QLS with a Mg + logic ion proceeds in the same way 
as with Be + , but the ground-state-cooling process [14[ 
has been modified to require only two lasers rather than 
three 15(. This cools the out-of- phase axial motional 
mode to an average quantum number of n < 0.05. It 



also enables quantum-non-demolition transfer (QNDT) 
of the Al + -clock-ion state to Mg + with approximately 
80 % fidelity in a single QNDT repetition and over 99 % 
fidelity after typically five QNDT repetitions fill ]. 

The trap utilized in the new A1+ clock has blade- 
shaped gold-coated beryllium-copper electrodes (Fig. Q} 
whose edges are approximately 400 [ira from the ions. 
The 25 Mg + logic ion (/ = 5/2) is manipulated with light 
of 279.5 nm wavelength from two independent lasers. A 
frequency-doubled dye laser resonantly drives | 2 S'i/2, F G 
{2, 3}) — ► 2 -p3/2 cycling and repumping transitions, while 
a frequency-quadrupled fiber laser is detuned by 40-60 
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TABLE I: Systematic effects that shift the clock from its ideal 
unperturbed frequency. Shifts and uncertainties given are in 
fractional frequency units (Av/v). See text for discussion. 



Effect 


Shift Uncertainty 
(lCT 18 ) (1CT 18 ) 


Excess micromotion 


n 

-y 





Secular motion 


-16.3 


5 


Blackbody radiation shift 


-9 


3 


Cooling laser Stark shift 


-3.6 


1.5 


Quad. Zeeman shift 


-1079.9 


0.7 


Linear Doppler shift 





0.3 


Clock laser Stark shift 





0.2 


Background gas collisions 





0.5 


AOM freq. error 





0.2 


Total 


-1117.8 


8.6 



TABLE II: Motional modes of the Al-Mg ion pair. For each 
of the six normal modes, the oscillation frequency and zero- 
point motional amplitude for each ion is listed. The trap 
axis corresponds to z, and x, y are two orthogonal radial di- 
rections whose orientation is determined by the trap geome- 
try. Also shown are the calculated {fie) and measured (um) 
Doppler-cooled average motional quantum numbers, and the 
time-dilation (TD) per motional quantum as well as the total 
TD per mode. 



f [MHz] 


6.53 


5.66 


5.20 


4.64 


3.41 


3.00 


25 Mg+ [am] 


4.9y 


2.9y 


4.65 


5.5i; 


-4.2i 


5.6z 


27 A1+ [nm] 


2.6y 


-5.00 


-4.12 


3.5i 


6.2x 


5.8z 


n M 


2.9 


4.5 


3.4 


6.3 


10.0 


7.0 


n c 


3.3 


3.8 


3.4 


5.9 


8.0 


5.9 


TD/quantum [10~ 18 ] 


0.226 


0.731 


0.197 


0.290 


0.771 


0.133 


Total TD [10 -18 ] 


0.77 


3.66 


0.77 


1.97 


8.10 


1.00 



GHz to drive Raman transitions between the 25 Mg + hy- 
pcrfine qubit states. Both ions are created via multi- 
photon-ionization of neutral atoms from ovens. A 396 
nm diode laser produces 27 A1 + 17J, while a frequency- 
doubled dye laser at 285 nm creates 25 Mg+ Typi- 
cally the clock transition is probed with 150 ms duration 
7r- pulses with a duty cycle of approximately 65 %. The 
remaining 35 % is occupied with state preparation and 
state detection functions as well as interleaved experi- 
ments that allow real-time measurement and reduction 
of micromotion. 

Systematic shifts of the Al-Mg clock are listed in 
Table |TJ Two types of residual motion cause time- 
dilation shifts: micromotion near the trap drive fre- 
quency of VRF = 59 MHz, and harmonic-oscillator (secu- 
lar) motion at the ion's normal mode frequencies (Ta- 
ble HI)) . In both cases the clock frequency shifts by 

¥ = -& (* + (iwmrJ 2 )' whcrc wc add to the rel " 

ativistic time-dilation (v 2 )/(2c 2 ) a frequency-dependent 
term that corresponds to the Stark shift from the motion- 
inducing electric fields. Here (v 2 ) is the mean-squared 
Al + -ion velocity and / is the frequency of motion. For 
the highest motional frequencies the Stark shift correc- 
tion is 2 %. 

Excess micromotion (EMM) refers to the rapid ion mo- 
tion at VRF It is caused by electric fields that force 
the ion away from the RF-minimum of the ion trap, 
or phase shifts between trap electrodes that cause the 
RF-fields to be non-zero at the pseudo-potential min- 
imum. We measure the amplitude of this motion at 
v rf by observing the motional-sideband strength of the 
Al + 1 So — > 3 Pi transition in three orthogonal direc- 
tions. For small amplitude of motion, the time dila- 
tion shift is Avjv = -\wrf/vl\ 2 = -2.8 x 10 15 1 | 2 , 
where Vl = 1.12 PHz is the probe laser frequency, and 
V = (^l i ^2,^3) is the measured EMM Lamb-Dickc pa- 
rameter (the ratio of sideband and carrier Rabi rates) 



for the three orthogonal directions. Typical values for 
Vi,2,3 are 0.01 to 0.04. We note that this EMM measure- 
ment method detects slow electric field fluctuations such 
as those caused by the migration of photo-electrons, as 
well as faster fluctuations caused by periodic line noise 
(50 to 60 Hz). Fluctuations that are shorter than the 
laser probe period of 0.05 to 0.1 ms will not be detected. 

In a perfect linear Paul trap EMM along the trap axis 
does not occur, but imperfections in the ion-trap geom- 
etry can lead to axial EMM . For the Al-Mg trap there 
exists a sharp minimum of axial EMM at one spatial 
location. The Al + ion is held at this point, but ran- 
dom background-gas collisions cause the Al-Mg ion pair 
to spontaneously re-order, which causes the Al + ion to 
move by 3 /im every 200 s, on average. The Al + ion is 
maintained at the position of minimal micromotion by 
adjusting the electrode voltages every 10 s to force the 
ion-pair into the desired order. When the ions are in the 
wrong order (about 5 % of the time), the Al + ion ex- 
periences excess axial micromotion and a corresponding 
clock shift of Av/v = -2.7 x 10~ 17 . This additional shift 
is included in Table HI 

During each Al + clock interrogation pulse, the Mg + 
ion is simultaneously Doppler cooled by a laser that is 
tuned 21 MHz below the \ 2 S 1/2 ,F = 3, to = -3) ->■ 
\ 2 Ps/2iF = 4, to = —4) cycling transition. The ampli- 
tude of secular motion (corresponding to the motional 
temperature) is extracted from the ratio of amplitudes for 
the red- and blue-sidebands of a Raman transition for all 
six normal modes [13] • Measured values are shown in Ta- 
ble [Hi together with values calculated from laser-cooling 
theory for a single ion [2(J. The single- Mg + -ion Doppler 
cooling limit is also valid when applied to each of the 
six normal modes of the Al-Mg ion pair. The measured 
and calculated motional quantum numbers agree within 
the measurement uncertainty, and we consider the stated 
30 % uncertainty for this shift to be a conservative limit. 

The Mg + Doppler cooling laser beam maintains the 
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Al-Mg ion pair at a constant motional temperature dur- 
ing the clock interrogation pulse, but because it also 
overlaps the Al + ion, it causes an AC Stark shift by 
coupling off-rcsonantly to allowed transitions that con- 
nect to the ground (3s 2 ) 1 So and excited (3s3p) 3 Po clock 
states. Following the evaluation of the blackbody radi- 
ation shift [Hj ], we estimate the differential clock polar- 
izability at 279.5 nm as Avjv = (-3.5 ± 0.6) x 10~ 17 S, 
where (S = I /Is) is the saturation parameter for 25 Mg+ 
(Is « 2470 W/m 2 ). The intensity I of the Mg+ Doppler 
cooling laser is estimated from the rate at which this 
laser repumps the | 2 5*1/2, F — 2,m = —2) dark hyper- 
fine ground state. The ion fluorescence photo-multiplier 
counts F(t) collected in a duration t due to repumping of 
the dark state may be written as F(t) = b(t+T(e~ t / T — 1)) 
where r = (0.217/S) ms and b is the bright-state count- 
ing rate. We extract S by fitting F(t) to the observed 
ion fluorescence. Typically we measure t = 2.1 ±0.8 ms, 
and find Avjv = (-3.6 ± 1.5) x 10~ 18 . 

Another Stark shift is caused by thermal blackbody ra- 
diation. The temperature of the Al-Mg ion trap is mea- 
sured with two platinum sensors at opposite ends of the 
trap structure which are expected to be at temperature 
extremes. Heat is removed primarily through thermal 
conduction at one end of the trap where we measure a 
temperature of 35 °C. At the higher thermal resistance 
trap-end the temperature is 40 °C, and the thermal ra- 
diation field impinging upon the ion is bounded by this 
maximum temperature and the laboratory room temper- 
ature of 22 °C: T lon = (31 ± 9) °C. 

During operation, the Al + clock servo alternates be- 
tween probing of the I 1 So, rap = 5/2)<H>| 3 P , m F = 5/2) 
and I 1 So, "IF = — 5/2)<-s>| 3 P , m F = — 5/2) transitions 
every 5 s, and the apparatus synthesizes an average of 
these two frequencies to eliminate first-order Zeeman 
shifts [il [2l| . Each transition's resonance is probed sev- 
eral times at the high- and low-frequency half-maximum 
points to derive a frequency-correction signal. The 
frequency-difference between the transitions is propor- 
tional to the mean magnetic field (£?), which allows an 
accurate estimate of the quadratic Zeeman shift due to 
the quasi-static quantization field of typically (B) =0.1 
mT. However, the quadratic Zeeman shift is proportional 
to (B 2 ) = (B) 2 + B\ c , where B\ G is the variance of the 
magnetic field about its mean. The dominant sources of 
varying magnetic fields are currents at vhf in conductors 
near the ion. We vary the trap RF drive power P and 
measure the frequency of the hyperfine clock transition 
in 25 Mg+ \F = 3,m F = 0) -> \F = 2,m F = 0) near 
1.789 GHz, which has a strong quadratic dependence on 
the magnetic field, to find B 2 AC = 1.45 x 10^ 12 (P/W) T 2 . 
For A1+ clock operation P = 15 W, and B\ G = 2.17 x 
10~ n T 2 , which alters the quadratic Zeeman shift by 
Avjv = (-1.4 ±0.3) x 10~ 18 . 

Other potential systematic shifts are listed in Table fl] 
When stabilizing the clock laser to the ion we probe 
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FIG. 2: Clock stability. Fractional frequency uncertainty vs. 
averaging period (r) for a comparison between the two Al + 
clocks (10700 s duration). Overlapping Allan deviation and 
N-sample standard deviation are shown [22] . For each com- 
parison measurement the coefficient of the r _1//2 asymptote 
is estimated and used to derive the measurement's statistical 
uncertainty. The 2.8 x 10~ 15 t~ 1//2 asymptote is reached for 
averaging periods that are longer than the servo time constant 
of 10 s. 

the clock transition alternately with counter-propagating 
laser beams to observe and cancel potential first-order 
Doppler shifts Q. For the clock comparison described 
below, we observe a differential shift for the two probe 
directions of (1.2 ± 0.7) x 10~ 17 but this effect is sup- 
pressed by taking the average. The suppression factor 
is limited because the atomic line-shape and hence the 
servo gains differ slightly for the two probe directions. 
During clock operation the mean fractional gain imbal- 
ance was 1.5 %, thereby reducing the possible first order 
Doppler shift to 3 x 10~ 19 . 

We have looked for Stark shifts due to the clock pulse 
itself by raising the intensity of the probe beam in one 
clock, and comparing the frequency to that produced by 
the other Al + clock. With an increase of 40 dB in the 
clock probe intensity, we observe no statistically signif- 
icant frequency difference with a fractional uncertainty 
of 2 x 10~ 15 , constraining any effect that scales linearly 
with this intensity to 2 x 10 -19 . We observe a rate of 
background gas collisions similar to the Al-Bc clock and 
assign a 5 x 10 -19 uncertainty to this potential shift 
Thermally induced frequency errors from beam switching 
AOMs (Fig. [TJ have been evaluated previously [Ej], and 
this uncertainty is reduced to 2 x 10~ 19 by operating the 
AOMs at less than 1 mW. 

We directly compared the Al-Mg clock with the Al- 
Be clock to perform independent tests that could reveal 
unaccounted-for clock shifts. The Al-Be clock was eval- 
uated with an accuracy of2.3xl0 -17 Q, and this eval- 
uation remained valid during the two-clock comparison. 
The two clocks were compared in 56 separate measure- 
ments, each of duration 1000 — 11000 s. As shown in 
Fig. [TJ the frequency of the 1 So0 3 P probe laser was 
actively steered to the Al + ion in the Al-Mg apparatus 
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FIG. 3: Measurements of the fractional frequency difference 
between the two Al + clocks (blue points). Error bars repre- 
sent the statistical uncertainty (see Fig. [2]). The horizontal 
line shows the weighted mean of —1.8 x 1CF 17 with an overall 
statistical uncertainty of ±7.0 x 1CP 18 (shaded band). 



with a servo time-constant of about 10 s. A portion of 
this laser light simultaneously probed the Al + ion in the 
Al-Bc clock, where it was servoed to the clock transi- 
tion in a separate digital feedback loop. The frequency 
produced by this feedback loop represents the frequency 
difference between the two clocks, and was recorded and 
analyzed for stability (Fig. [5]). Average frequencies of 
the individual measurements corrected for known shifts 
are shown in Fig. [3J where the overall weighted mean is 
{vALMg - v MBe )/v = (-1.8 ± 0.7) x 10" 17 . This value is 
consistent with the 1-er error of 2.5 x 10 -17 that is calcu- 
lated by adding in quadrature the inaccuracies of the two 
clocks and the statistical uncertainty. The reduced-^ 2 for 
this data set is 1.02, indicating that the error bars derived 
from the estimated stability correctly capture the scat- 
ter of the data. For this data set the total in-loop servo 
error for the Al-Mg clock was Av/v = 6 x 10 -19 , and 
for a well-designed servo loop this error declines faster 
than the statistical error. It is therefore not included as 
a systematic shift in Table fl] 

In summary, we have built an Al + ion clock with a 
fractional frequency inaccuracy of 8.6 x 10 -18 . Its fre- 
quency is compared to that of a previously constructed 
Al + clock, and the measured fractional frequency differ- 
ence of (— 1.8 ± 0.7) x 10~ 17 is consistent with the inac- 
curacy of the previous clock. Significantly, the statistical 
uncertainty in the frequency comparison of 7.0 x 10~ 18 is 
smaller than the inaccuracy of either clock, and the av- 
erage measurement stability was 2.8 x 10 _15 r -1 / 2 (total 



measurement duration: 164 967 s). This result may be 
compared to other direct same-species atomic clock com- 
parisons, where two Cs microwave clocks have reached 
an agreement of (4 ±3) x 10 -16 [23j, and two Yb + single 
ion clocks showed agreement of (3.8 ± 6.1) x 10 -16 [24|. 



Future frequency ratio measurements of the Al + clock 
and the NIST Hg + optical clock would enable improved 
constraints on present-era changes in the fine-structure 
constant [f|. 
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